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Abstract — Charge collection measurements performed on heav- 
ily irradiated p-spray DOFZ pixel sensors with a grazing angle 
hadron beam provide a sensitive determination of the electric field 
within the detectors. The data are compared with a complete 
charge transport simulation of the sensor which includes free 
carrier trapping and charge induction effects. A linearly varying 
electric field based upon the standard picture of a constant type- 
inverted effective doping density is inconsistent with the data. A 
two-trap double junction model implemented in the ISE TCAD 
software can be tuned to produce a double peak electric field 
which describes the data reasonably well. The modeled field differs 
somewhat from previous determinations based upon the transient 
current technique. The model can also account for the level of 
charge trapping observed in the data. 

Index Terms — Pixels; Radiation effects; Space charge; Simula- 
tion; Electric fields; 



I. Introduction 

IN the recent years, detectors consisting of one and two 
dimensional arrays of silicon diodes have come into 
widespread use as tracking detectors in particle and nuclear 
physics experiments. It is well understood that the intra-diode 
electric fields in these detectors vary linearly in depth reaching 
a maximum value at the p-n junction. The linear behavior 
is a consequence of a constant space charge density, N e g, 
caused by thermodynamic ally ionized impurities in the bulk 
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material. It is well known that the detector characteristics are 
affected by radiation exposure, but it is generally assumed 
that the same picture is valid after irradiation. In fact, it is 
common to characterize the effects of irradiation in terms of a 
varying effective charge density. The use of N e g to characterize 
radiation damage has persisted despite a growing body of 
evidence [1]- [6] that the electric field does not vary linearly as 
a function of depth after heavy irradiation but instead exhibits 
maxima at both n + and p + implants. The study presented in 
this paper demonstrates conclusively that the standard picture 
does not provide a good description of irradiated silicon pixel 
detectors. We show that it is possible to adequately describe 
the charge collection characteristics of a heavily irradiated 
silicon detector in terms of a tuned double junction model 
which produces a double peak electric field profile across the 
detector. The allowed parameter space of the model can also 
accommodate the expected level of leakage current and the level 
of charge trapping observed in the detector. 

This paper is organized as follows: Section UTI describes the 
experimental technique and data, Section [HI] describes the car- 
rier transport simulation used to interpret the data, Section llVl 
describes the technique used to model double peak electric 
fields and the limitations of previous models. The tuning of 
a successful model is discussed in Section [V] Section IV1I 
summarizes the results and develops several conclusions. 

II. Experimental Technique and Data 



This investigation is based upon beam test data that were 
accumulated as part of a program to develop a silicon pixel 
tracking detector [7] for the CMS experiment at the CERN 
Large Hadron Collider. The measurements were performed 
in the H2 line of the CERN SPS in 2003/04 using 150-225 
GeV pions. The beam test apparatus is described in [8] and 
is shown in Fig. ^ A silicon beam telescope [9] consisted of 
four modules each containing two 300 /mi thick single-sided 
silicon detectors with a strip pitch of 25 /im and readout pitch 
of 50 /im. The two detectors in each module were oriented to 
measure horizontal and vertical impact coordinates. A pixel 
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hybrid detector was mounted between the second and third 
telescope modules on a cooled rotating stage. A trigger signal 
was generated by a silicon PIN diode. The analog signals from 
all detectors were digitized in a VME-based readout system 
by two CAEN (V550) and one custom built flash ADCs. The 
entire assembly was located in an open-geometry 3T Helmholtz 
magnet that produced a magnetic field parallel or orthogonal to 
the beam. The temperature of the test sensors was controlled 
with a Peltier cooler that was capable of operating down to 
-30°C. The telescope information was used to reconstruct the 
trajectories of individual beam particles and to achieve a precise 
determination of the particle hit position in the pixel detector. 
The resulting intrinsic resolution of the beam telescope was 
about 1 /im. 

Beam telescope modules 
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Fig. 1 

Side view of the beam test apparatus consisting of four 
horizontal and four vertical planes of silicon strip detectors 
and a rotating stage for the pixel detector. 



after the deposition of under-bump metalization and indium 
bumps. A number of sensors were irradiated at the CERN PS 
with 24 GeV protons. The irradiation was performed without 
cooling or bias. The delivered proton fluences scaled to 1 MeV 
neutrons by the hardness factor 0.62 [11] were 6x 10 14 n cq /cm 2 
and 8xl0 14 n eq /cm 2 . The former sample was annealed for 
three days at 30°C. In order to avoid reverse annealing, the 
sensors were stored at -20°C after irradiation and kept at 
room temperature only for transport and bump bonding. All 
sensors were bump bonded to PSI30/AC30 readout chips [12] 
which allow analog readout of all 704 pixel cells without zero 
suppression. The PSI30 also has a linear response to input 
signals ranging from zero to more than 30,000 electrons. 

B. Data 

The main focus of the work presented in this paper involves 
a set of charge collection measurements that were performed 
using the "grazing angle technique" [13]. As is shown in Fig. 12 
the surface of the test sensor is oriented by a small angle 
(15°) with respect to the hadron beam. A large sample of 
data is collected with zero magnetic field and at a temperature 
of — 10°C. The charge measured by each pixel along the y 
direction samples a different depth z in the sensor. Precise entry 
point information from the beam telescope is used to produce 
finely binned charge collection profiles. For unirradiated sen- 
sors, the cluster length determines the depth over which charge 
is collected in the sensor. 
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A. Pixel Hybrids 

The prototype pixel sensors are so-called "n-in-n" devices: 
they are designed to collect charge from n + structures im- 
planted into n- bulk silicon. This design is thought to provide 
greater longevity in high radiation fields and to allow "partially- 
depleted" operation after "type-inversion" of the substrate. It 
is more expensive than the "p-in-n" process commonly used 
in strip detectors because it requires double-sided processing 
and the implementation of inter-pixel isolation. Two isolation 
techniques were tested: p-spray, where a uniform medium dose 
of p-impurities covers the whole structured surface, and p-stop, 
where higher dose rings individually surround the n + implants. 
Results on the Lorentz angle and charge collection efficiency 
measurements as well as a detailed description of both designs 
can be found elsewhere [8], [10]. In this paper we discuss only 
measurements performed on p-spray sensors. All test devices 
were 22x32 arrays of 125x125 /im 2 pixels having a sensitive 
area of 2.75x4 mm 2 . The substrate was 285 /im thick, n-doped, 
diffusively-oxygenated silicon of orientation (111), resistivity 
2-5 k£!-cm and oxygen concentration in the order of 10 17 cm -3 . 
Individual sensors were diced from fully processed wafers 



Fig. 2 

The grazing angle technique for determining charge 
collection profiles. the cluster length is proportional to the 
depth over which charge is collected. 



The profiles that were observed for an unirradiated sensor 
and for a sensor that was irradiated to a fluence of $ = 
8 x 10 14 n eq /cm 2 are shown in Fig. [3] as function of the 
distance from the beam entry point. The unirradiated sensor 
was operated at a bias voltage of 150 V which is well above 
its depletion voltage (approximately 70 V). It produces a large 
and uniform collected charge distribution indicating that it is 
fully depleted (a large y coordinate indicates a large collection 
distance). The irradiated sensor was operated at bias voltages 
varying from 150 V to 600 V. It appears to be partly depleted at 
150 V, however, signal is collected across the entire thickness 
of the detector. Another puzzle is that the ratio of the charges 
collected at 300 V bias and 150 V bias and integrated along 
the distance to the beamentry point is 2.1 which is much larger 
than the maximum value of \[2 expected for a partially depleted 



junction, where the depletion depth is proportional to the square 
root of the bias voltage [14]. It is clear that the profiles for 
the irradiated sensor exhibit rather different behavior than one 
would expect for a heavily-doped, unirradiated sensor. 
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Charge collection profiles for an irradiated 

(<J) = 8 X 10 14 N cq /cm 2 ) AND AN UNIRRADIATED SENSOR 
($ = Neq/cm 2 ) OPERATED AT SEVERAL BIAS VOLTAGES. 



III. Simulation and comparison with data 

It is well-known that carrier trapping is a significant effect 
in heavily irradiated silicon detectors. In order to evaluate the 
effects of trapping, it is necessary to implement a detailed 
simulation of the sensor. Our simulation, PIXELAV [15], [16], 
incorporates the following elements: an accurate model of 
charge deposition by primary hadronic tracks (in particular to 
model delta rays); a realistic electric field map resulting from 
the simultaneous solution of Poisson's Equation, continuity 
equations, and various charge transport models; an established 
model of charge drift physics including mobilities, Hall Effect, 
and 3-D diffusion; a simulation of charge trapping and the sig- 
nal induced from trapped charge; and a simulation of electronic 
noise, response, and threshold effects. A final step reformats the 
simulated data into test beam format so that it can be processed 
by the test beam analysis software. 

Several of the PIXELAV details described in [15], [16] 
have changed since they were published. The commercial 
semiconductor simulation code now used to generate a full three 
dimensional electric field map is the ISE TCAD package [17]. 
The charge transport simulation was modified to integrate only 
the fully-saturated drift velocity, 



dr 
dt 



\x qE + fiR H E x B + qfi 2 R 2 H (E ■ B)B 



1 + fi 2 R 2 H \B 



(1) 



where n(E) is the mobility, q = ±1 is the sign of the charge 
carrier, E is the electric field, B is the magnetic field, and 
Rh is the Hall factor of the carrier. The use of the fully- 
saturated drift velocity permits much larger integration steps 



(which had previously been limited by stability considerations) 
and significantly increases the speed of the code. 

The simulation was checked by comparing simulated data 
with measured data from an unirradiated sensor. A plot of the 
charge measured in a single pixel as a function of the horizontal 
and vertical track impact position for normally incident tracks is 
shown in Fig.|4] The simulation is shown as the solid histogram 
and the test beam data are shown as solid points. Note that 
the sensor simulation does not include the "punch-through" 
structure on the n + implants which is used to provide a high 
resistance connection to ground and to provide the possibility 
of on-wafer IV measurements [18]. There is reduced charge 
collection from this portion of the implant and the data show 
reduced signal in both projections at the bias dot. Another 
check, shown in Table HJ is the comparison of the average 
Lorentz angle measured at several bias voltages [8]. In both 
cases, reasonable agreement is observed between measured and 
simulated data. 
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Fig. 4 



Collected charge measured in a single pixel as a function of 
the horizontal (a) and vertical (b) track impact position for 
tracks that are normally incident on an unirradiated sensor, 
the simulation is shown as a solid histogram and the test beam 
data are shown as solid dots. 



TABLE I 

Measured and simulated values of average Lorentz angle Q l 
versus bias voltage for an unirradiated sensor 



Bias Voltage [V] 


Measured 8l [deg] 


Simulated 8l [deg] 


150 


22.8 ±0.7 


24.7±0.9 


300 


14.7 ±0.5 


17.4±0.9 


450 


11.2 ±0.5 


12.0±0.9 



The charge collection profiles for a sensor irradiated to a 
fluence of $ = 6 x 10 14 n eq /cm 2 and operated at bias voltages 



of 150 V and 300 V are presented in Fig. [5(a)| and |5(Fj| 
respectively. The measured profiles are shown as solid dots and 
the simulated profiles are shown as histograms. The simulated 
profiles were generated with electric field maps corresponding 
to two different effective densities of acceptor impurities. The 
full histograms are the simulated profile for N c s — —4.5 x 



10 cm . Note that the 300 V simulation reasonably agrees 
with the measured profile but the 150 V simulation is far too 
broad. The dashed histograms show the result of decreasing 
iV e ff to —24 x 10 12 cm~ 3 . At this effective doping density, the 
width of the simulated peak in the 150V distribution is close to 
correct but it does not reproduce the second peak observed in 
the data at large y. The 300 V simulated distribution is far too 
narrow and the predicted charge is lower than the data (note 
that the profiles are absolutely normalized). It is clear that a 
simulation based upon the standard picture of a constant density 
of ionized acceptor impurities cannot reproduce the measured 
profiles. 
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The measured and simulated charge collection profiles for a 
sensor irradiated to a fluence of $ = 6 x 10 14 n oq /cm 2 . the 

profiles measured at bias voltages of 150 v (a) and 300 v (b) are 
shown as solid dots. the full histograms are the simulated 

profiles for a constant effective doping a^ff 



-4.5 X 10 12 CM 



OF ACCEPTOR IMPURITIES. THE DASHED HISTOGRAMS ARE THE 
SIMULATED PROFILES FOR A CONSTANT EFFECTIVE DOPING 
7V cff = -24 X 10 12 CM" 3 . 



Note that the simulation of this irradiated sensor includes the 
effects of trapping. The trapping rates of electron and holes have 
been shown to scale linearly with fluence [19]. Unfortunately, 
the measured fluences have a fractional uncertainty of ±10% 
which feeds directly into the uncertainty on the trapping rates. 
Additional uncertainty arises because annealing can modify the 
trapping rates by 30% [19] leading to a fairly large overall 
uncertainty. 

IV. Two-trap Models 

The large number of measurements suggesting that large 
electric fields exist at both sides of an irradiated silicon diode 
has given rise to several attempts to model the effect [3], [6], 
[20]. The most recent of these by Eremin, Verbitskaya, and Li 
(EVL) [20] is based upon a modification of the Shockley-Read- 
Hall (SRH) statistics. The EVL model produces an effective 



space charge density p e g from the trapping of free carriers 
originated from the leakage current by one acceptor trap and 
one donor trap. The effective charge density is related to the 
occupancies and densities of traps as follows, 



Peg = e \N D f D - Na/a] + Pdopants 



(2) 



where: Np and Na are the densities of donor and acceptor trap- 
ping states, respectively; fr> and Ja are the occupied fractions 
of the donor and acceptor states, respectively, and Pdopants is 
the charge density due to ionized dopants. Charge flows to and 
from the trapping states due to generation and recombination. 
The occupied fractions are given by the following standard SRH 
expressions, 

a E D /kT 

d[3) 



fn 

fA 

where: 
holes, 



v e a®n,.<: ' 



v h <T A nie~ 



E A /kT 



v e a A (n + me EA / kT ) + v h a A (p + me- EA l kT ) 
v e and Vh are the thermal velocities of electrons and 



respectively; a® and ajf are the electron and hole 



and a A 



are 



capture cross sections for the donor trap; a J e 
the electron and hole capture cross sections for the acceptor 
trap; n and p are the densities of free electrons and holes, 
respectively; rii is the intrinsic density of carriers; Ed and 
Ea are the activation energies (relative to the mid-gap energy) 
of the donor and acceptor states, respectively. The generation- 
recombination current caused by the SRH statistics for single 
donor and acceptor states is given by the following expression, 



U 



v h v e a%a^N D (np - nf) 



'(n + me ED / kT ) + v h a®(p + me- ED / kT ) 



v h v e a A a A N A (np- 



n?) 



l (n + me E */ kT ) + v h a A (p + n i e- £ W cT ) 



•(4) 



Within the EVL model, the four trapping cross sections are 
set to 10~ 15 cm 2 . The leakage current is generated from an 
additional SRH trapping state that is introduced for this purpose 
but is assumed not to trap charge. The donor and acceptor 
states are assumed not to generate leakage current which, 
given the small size of the cross sections, is a self-consistent 
assumption. The densities of the donor and acceptor states (Np 
and Na) are adjusted to "fit" data obtained with the Transient 
Current Technique (TCT). The parameters of the model are 
given in Table ITT1 The trap densities are scaled to fluence and 
are given in terms of introduction rates g mt = NA/o/^cq- 
An illustrative sketch of the EVL model has been reproduced 

TABLE II 

Parameters of the EVL Model [20]. E v and Eq are the valence 

AND CONDUCTION BAND ENERGY LEVELS, RESPECTIVELY 



Trap 


E (eV) 


flint (cm -1 ) 


a e (cm z ) 


cr h (cm 2 ) 


Donor 


E v + 0.48 


6 


1 X 10" ib 


1 x 10" ib 


Acceptor 


E c - 0.525 


3.7 


1 X 10~ 15 


1 x 10~ 15 



from [20] and is shown in Fig. [6] Figure |6ja) shows a uniform 



current density flowing across a reverse-biased junction. Since 
holes are produced uniformly across the junction and flow to 
the p + backplane, the hole current density increases linearly 
with increasing z from the n + implant to the p + implant. 
The electrons flow to the n + implant and the electron current 
density increases with decreasing z. The actual carrier densities 
depend upon the details of the fields and mobilities but vary 
monotonically across the junctions as shown in Fig. EJb). The 
trapping of the mobile carriers produces a net positive space 
charge density near the p + backplane and a net negative space 
charge density near the n + implant as shown in Fig. |6jc). 
Since positive space charge corresponds to n-type doping and 
negative space charge corresponds to p-type doping, there 
are p-n junctions at both sides of the detector. The electric 
field in the sensor follows from a simultaneous solution of 
Poisson's equation and the continuity equations. The resulting 
z-component of the electric field is shown in Fig.|6jd). It varies 
with an approximately quadratic dependence upon z having a 
minimum at the zero of the space charge density and maxima 
at both implants. 




Fig. 6 

An illustrative sketch of the EVL model for a reverse biased 

DEVICE [20]. 



In order to test whether the electric field predicted by 
the EVL model would improve the agreement between the 
simulated and measured charge collection profiles shown in 
Section [H]] it was necessary to implement the EVL model 
in ISE TCAD. TCAD contains a complete implementation of 
SRH statistics. However, the EVL modifications of SRH are 
not incorporated. In particular, any state added to generate 
leakage current would also trap charge. It is possible to use the 
TCAD "Physical Model" interface to replace the existing SRH 
implementation with a modified one, however, a less invasive 
approach was adopted. 

Our approach is based upon three observations: 

1) The trapping cross sections are poorly known. The cross 
sections for states observed in various types of defect 
spectroscopy vary over several orders of magnitude. 

2) The occupancies fo/A of the trapping states are inde- 
pendent of the scale of the cross sections. If the capture 



cross sections cr e /h in Eq. [5] are rescaled by a factor r, 
then f D i& are unchanged. The occupancies depend only 
upon the ratio of electron and hole cross sections, ah/a e , 
3) The current U generated by the donor and acceptor 
impurities is linear in the cross section rescaling factor 
r. 

These observations imply that it is possible to implement the 
EVL model in TCAD simply by setting erf = erf = = afi 
and by varying the size of the common cross section until the 
generation current is equal to the observed or expected leakage 
current. The trap occupancies are not affected in zeroth order 
by the rescaling, but the leakage current and the free carrier 
densities are affected by r. The carrier densities have a first- 
order effect on the occupancies so that varying r does alter p c s ■ 
This approach uses the same trapping states to produce space 
charge and leakage current (it is not necessary to introduce 
current-generating states). 

Using the activation energies and introduction rates for the 
donor and acceptor states given in Table |ll| the only free 
parameter in the TCAD implementation of the EVL model 
is the size of the common cross section or equivalently, the 
leakage current. The leakage current in the test sensors was 
observed to increase substantially after bump-bonding to the 
readout chips. Although the cause of the increased current is 
not clear, it is possible that it is due to increased surface/edge 
leakage or that it is caused by stressing the detector. The 
leakage current of the sample irradiated to 6 x 10 14 n eq /cm 2 
at a bias voltage of 300 V and temperature of -10°C was 15.1 
/j, A. The expected value calculated using Eq. 1 and 2 of [21] 
and a damage rate ao = 4 x 10~ 17 A/cm is found to be 4.7 
fiA. 

The EVL model was implemented in the sensor simulation 
using several parameter choices. The resulting charge collection 
profiles for a sensor irradiated to a fluence of $ = 6 x 
10 14 n cq /cm 2 and operated at 150 V, 200 V and 300 V are 
shown in Fig. The measured profiles are again shown as 
solid dots. The solid histogram shows the EVL model with the 
leakage current adjusted to the expected value by setting <J e /h = 
0.47 x 10 -14 cm 2 . It clearly substantially underestimates the 
collected charge at both voltages. The effect of increasing 
the leakage current to the measured value is shown as the 
dashed histogram, where the trapping cross sections are set 
to a e /h = 1.48 x 10~ 14 cm 2 . This increases the signal at 300 
V but does not reproduce the shape of the high z tail of the 
150 V distribution. We conclude that the EVL model does not 
describe the measured charge collection profiles. 

V. An Improved Two-Trap Model 

Although the EVL model does not describe the measured 
charge collection profiles, the electric field generated by the 
two-trap mechanism could potentially reproduce the main fea- 
tures observed in the data. At low bias voltages, the combi- 
nation of the quadratic minimum in the electric field and free 
carrier trapping can act like a "gate" suppressing the collection 
of charge from the p + side of the detector. The measured profile 
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Fig. 7 

THE MEASURED CHARGE COLLECTION PROFILES FOR A SENSOR IRRADIATED TO A FLUENCE OF <E> = 6 X 10 14 N oq /CM 2 (SOLID DOTS) AND OPERATED AT 

150 V (A), 200 V (B) AND 300 V (C) ARE COMPARED WITH SIMULATIONS BASED UPON THE EVL MODEL. THE SOLID HISTOGRAM SHOWS THE EVL 
MODEL WITH THE LEAKAGE CURRENT ADIUSTED TO THE EXPECTED VALUE. THE DASHED HISTOGRAM SHOWS THE EFFECT OF INCREASING THE LEAKAGE 

CURRENT TO THE MEASURED VALUE. 



would then appear to be a narrow peak on the n + side of the 
detector. As the bias is increased, the magnitude of the field at 
the minimum would also increase and effectively "lift the gate" 
which would allow much more charge collection from the p + 
side of the detector. 

In order to investigate whether a two-trap double junction 
EVL-like model can describe the measured charge collection 
profiles, a tuning procedure was adopted. Relaxing the EVL 
requirement that all trapping cross sections are equal, the model 
has six free parameters (Nd, Na, cr^, cr^, erf, af) that 
can be adjusted. The activation energies are kept fixed to the 
EVL values. Additionally, as was discussed in Section [HI] the 
electron and hole trapping rates, T e and Th, are uncertain at the 
30% level due to the fluence uncertainty and possible annealing 
of the sensors. They are treated as constrained parameters. The 
parameters of the double junction model were systematically 
varied and the agreement between measured and simulated 
charge collection profiles was judged subjectively. 

In the course of the tuning procedure, it became clear that 
the EVL model does not produce a sufficiently large electric 
field on the p + side of the detector. The solution to this 
problem is to increase the density of donors (hole traps) as 
compared to the density of acceptors (electron traps). When 
this is done, the z position of the minimum in the effective 
charge density shifts toward the n + implant as sketched in 
Fig. EJa). Unfortunately, this causes the "peak" in the 150 V 
simulated charge profile to become too narrow. The position of 
the charge density minimum can be restored to a position nearer 
the midplane of the detector by decreasing the ratios of the hole 



and electron cross sections from 1.0 to 0.25 (o/^/cf = 0.25 
and erf /af = 0.25) as shown schematically in Fig.|8jb). Note 
that the adjustment of the cross sections for both trap types 
minimizes the field in the quadratic minimum while allowing 
for large fields at the implants. For simplicity it is assumed that 
the electron cross sections are equal (a^ = erf = a e ) and that 
the hole cross sections are equal (cr® — af = (Th). 




Fig. 8 

The effect of increasing N d /Na when (a) the electron and 

HOLE CROSS SECTIONS ARE EQUAL, AND WHEN (B) cr h /cr e = 0.25. 



The current best "fit" to the measured charge collection 
profiles is called 'BF' and reduces the ratio between the 
densities of acceptor and donor states, Na/Nd, from the EVL 
value of 0.62 to 0.40. The z-component of the simulated electric 
field, E z , is plotted as a function of z in Fig.|9]for bias voltages 
of 150 V and 300 V. The field profiles have minima near 
the midplane of the detector. Note that the minimum field at 
150 V bias appears to be very small but is still approximately 
400 V/cm. The electric field profiles resulting from a constant 



p-type doping of density N e g = —4.5 x 10 12 cm 3 are shown 
as dot-dashed and dotted curves for reference. 
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The ^-component of the simulated electric field resulting from 
the model bf is plotted as a function of 2. the field profiles 
for 150 v and 300 v are shown as as solid and dashed curves, 

respectively. the electric field profiles resulting from a 
constant p-type doping of density n e ff = —4.5 x 10 12 cm -3 are 
shown as dot-dashed and dotted curves for 150 v and 300 v, 
respectively. 



The measured charge collection profiles at bias voltages 
between 150 V and 450 V are compared with the BF simulation 
in Fig. [K)] for a fiuence of 6 x 10 14 n oq /cm 2 . The electron 
trapping rate for the lower fiuence is set to 93% of the nominal 
value. 

Although for high values of the bias voltage the simulation 
falls below the measured profile, it provides a reasonable de- 
scription of the measurements. Several features of the measured 
distributions are described well by the simulation. Note that 
both data and simulations show a distinctly negative signal 
near y — ^m. This can be understood as a consequence 
of hole trapping. Electrons deposited near the n + implant are 
collected with high efficiency whereas holes deposited near 
the implant must transit the entire detector thickness to reach 
the p + backplane. If the holes are collected, they produce no 
net signal on the n + side of the detector. However if the 
holes are trapped, then a negative signal is induced and is 
most visible in the y < region. Another feature is the 
"wiggle" in the 150 V profiles. The relative signal minimum 
near y = 700 /im corresponds to the E z minimum where both 
electrons and holes travel only short distances before trapping. 
This small separation induces only a small signal on the n + side 
of the detector. At larger values of y, E z increases causing the 
electrons drift back into the minimum where they are likely to 
trap. However, the holes drift into the higher field region near 
the p + implant and are more likely to be collected. The net 
induced signal on the n + side of the detector therefore increases 
and creates the local maximum seen near y = 900 /im. 

The BF model fixes the ratios N A /No and the ratio between 



cross sections ah,ja e leaving the parameters No and a e to vary. 
Over a restricted range, an increase in No can be offset by a 
decrease in a e keeping the electric field profile approximately 
unchanged. Scaling the electron cross section as a e oc N^ 25 
produces very similar charge collection profiles. The allowed 
region in the No-cr e space is shown in Fig. II If a') as the solid 
line in the logarithmic space. If the donor density becomes too 
small (Njj < 20 x 10 14 cm -3 ), the 150 V simulation produces 
enhanced signal at large z. If the donor density becomes 
too large {No > 50 x 10 14 cm -3 ), the 300 V simulation 
produces insufficient signal at large z. Since the simulated 
leakage current varies as 7i oa k oc a e No, different points on the 
allowed solid contour correspond to different leakage current. 
Contours of constant expected leakage current are shown as 
dashed curves and are labeled in terms of the damage parameter 
a in units of a$. It is clear that the simulation can accommodate 
the expected leakage current which is smaller than the measured 
current by a factor of three. 
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THE ALLOWED REGION IN THE N D -a e SPACE FOR THE BF MODEL IS 
SHOWN AS THE SOLID LINE IN (A) AND (B). CONTOURS OF CONSTANT 
LEAKAGE CURRENT ARE SHOWN AS DASHED CURVES IN (A) AND ARE 
LABELED IN TERMS OF THE CORRESPONDING DAMAGE PARAMETER a IN 

units of ao. Contours of constant electron trapping rate are 

SHOWN AS DASHED CURVES IN (B) AND ARE LABELED IN TERMS OF THE 
UN-ANNEALED TRAPPING RATE Tq FOR THE NOMINAL FLUENCE. 



The electron and hole traps in the model should also con- 
tribute to the trapping of signal carriers. The contributions of 
these states to the effective trapping rates of electrons and holes 
are given by the following expressions 

T e = v e [afN A {\ - f A ) + a^Nofo] * v e afN A (5) 
T h = v h [a%N D {l - f D ) + e£N A f A ] ~ v h a£ N D (6) 

where it has been assumed that the trap occupancies are 
small and the thermal velocity of electrons at -10°C is set 
to v e = 2.15 x 10 7 cm/s. Because N A /No is assumed to be 
constant, contours of constant electron trapping rate are parallel 
to contours of constant leakage current in No-<J e space. The 
best "fit" of the simulation to the measured profiles reduced T e 
to 93% of the un-annealed trapping rate To — Q[3 e — 0.33 ns -1 
for the nominal fiuence $ [19]. These contours are compared 




with the allowed contour in Fig. II 1T b ), It is clear that the 
simulation can accommodate the measured trapping rate in the 
same region of parameter space that maximizes the leakage 
current. 

Figure II lf b) also suggests a solution to the puzzle that 
the trapping rates have been shown to be unaffected by the 
presence of oxygen in the detector bulk [19] whereas it is 
well-established that the space charge effects are quite sensitive 
to the presence of oxygen in the material [22], [23]. It is 
clear from Fig II li b) that small-cross-section trapping states 
can play a large role in the effective charge density but a 
small one in the effective trapping rates: every point on the BF 
line produces 100% of the effective charge density but only 
the larger cross section points contribute substantially to the 
trapping rate. If the formation of the additional small-cross- 
section states were suppressed by oxygen, then p e s could be 
sensitive to oxygenation whereas T e /h would be insensitive to 
oxygenation. This is another consequence of the observation 
that the occupancies f D ^ of the trapping states are independent 
of the scale of the cross sections in the steady state (see 
Section ITvV The trapping of free carriers is not a steady-state 
phenomenon and is sensitive to the scale of the trapping cross 
sections. 

VI. Conclusions 

The main result of the work presented in this paper is that 
a double peak electric field is necessary to describe the charge 
collection profiles measured in heavily irradiated pixel sensors. 
A simulation utilizing a linearly varying electric field based 
upon the standard picture of a constant type-inverted effective 
charge density is inconsistent with the measurements. 

A two-trap EVL-like model can be tuned to provide a 
reasonable description of the measurements. It can also account 
for the expected level of the leakage current (although not the 



observed leakage current) and the observed electron trapping 
rate. It is important to state that any two-trap model is, at 
best, an "effective theory". It is well-known that there are 
many trapping states in heavily irradiated silicon that trap 
charge. There may also be thermodynamically ionized defects 
that contribute to the effective space charge density. Clearly a 
two-trap model can describe the gross features of the physical 
processes in our sensors but it may not be able to describe all 
details. This also implies that the parameters of the two-trap 
model presented in this paper are unlikely to have physical 
reality. 

The charge-sharing behavior and resolution functions of 
many detectors are sensitive to the details of the internal 
electric field. A known response function is a key element of 
any reconstruction procedure. A working effective model will 
permit the detailed response of these detectors to be tracked as 
they are irradiated in the next generation of accelerators. 

Finally, we note that quantities like W (depletion depth) 
and iVcff, which are related to the picture of uniform type 
inversion in irradiated silicon sensors, may correctly suggest 
reduced detector performance but given the evidence of double 
peak electric fields and free carrier trapping, have no physical 
significance. 
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